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Abstract 

In an attempt to provide a clearer understanding of the impressive increase in Tc under pressure in elemental Li, linear response 
I calculation of the phonon dispersion curves, electron-phonon matrix elements, phonon linewidths and mode A's have been carried 
■ O out on a finer mesh (24^ in the Brillouin zone) than done previously, for the volume corresponding to 20 GPa pressure. The 
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result illustrates the great need for a fine mesh (even finer than this) for converged results of A and the spectral function a^F. 
Although the initial pressure-induced transverse Ti phonon instability (in harmonic approximation) near the symmetry point K 
jhas dominated attention, the current results show that the high value of Tc gets strong contributions from elsewhere in the zone, 
particularly from the longitudinal mode along (100). 
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1. Introduction 

■ Elemental lithium, with superconducting Tc ~ 15 K 
around[l,2,3] 30 GPa and up to 20 K in strained samples[l] 
at 50 GPa, joins MgBj (Tc=40 K), metallic Y (Tc=20 
K)[4], and PuCoGas (Tc=19 K)[5] as the big surprises in 
Critical temperature in this century. Transformation from 
a simple nearly- free-electron metal not superconducting [6] 
cibove 100 iiK to a strongly coupled electron-phonon super- 
conductor at 20-50 GPa has been explained, at least semi- 
quantitatively, in recent work by Profeta et al.[7] and by 
Kasinathan et al.[8] . The evolution of Tc has been linked to 
the increase in electron-phonon coupling of specific phonon 
branches along intersections of Kohn anomaly surfaces with 
pressure. The coupling strength is quite strongly dependent 
on the phonon wave vector Q and phonon polarization. In 
particular, the < 110 > polarized phonon along the (110) 
T-K direction softens with increasing pressure and finally 
becomes unstable. 

The nesting function S,{Q) which describes the phase 
space for the electron-hole scattering processes across the 
Fermi surface has fine structure [8] requiring a very dense 
mesh to obtain accurate values where its amplitude is 



large. Since the electron-phonon (EP) coupling involves 
these same scattering events weighted by EP matrix ele- 
ments that normally vary smoothly with Q, the fine struc- 
ture in ^(Q) implies corresponding fine structure in EP 
coupling (as contained for example in the EP contribution 
to the phonon linewidths 7qi/). Such fine structure makes 
the sampling of the Brillouin zone an issue if accurate 
numerical results for a^F are desired. Our previous linear 
response phonon calculations [8] were performed on 12^ 
mesh (72 Q points in the irreducible zone). In the present 
work, we focus on the 20 GPa case (fee lattice constant 
ap = 6.80 bohr) where strong phonon softening arises but 
structural instability has not yet occurred. The phonon en- 
ergies and the Eliashberg spectral function constants have 
been presently calculated on a much finer 24"^ mesh (413 Q 
points in the irreducible zone). Careful examination of the 
results lead to a revised picture of what coupling is impor- 
tant for Tc (as opposed to what is important in increasing 
A, which is a far simpler question). 



2. Computational Details 



* Corresponding author. 

Email address: pickettaphysics.ucdavis.edu (W. E. Pickett ). 



The electronic wave functions are calculated using a full 
potential linear muffin tin orbital (FP LMTO) method as 
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Fig. 1. (color online) Calculated phonon spectrum for fee Li at 20 
GPa along high symmetry directions. The solid diamonds denote 
the longitudinal mode £, and where the transverse branches are 
non-degenerate, solid circles denote the transverse mode 71 and open 
circles denote the transverse mode 72 . The mode symmetry has been 
used to connect branches across crossings (F-K) although this is not 
possible off symmetry lines where branches do not cross. 

implemented by Savrasov [10]. We have used the Vosko- 
Wilk-Nusair local exchange-correlation [11]. The phonon 
energies and the coupling constants have been calculated 
using linear response theory [10]. Phonons were calculated 
for a dense 24'^ mesh corresponding to 413 inequivalent 
Q points in the irreducible zone. The electronic BZ zone 
integration was performed with a 40^ k point grid. 



3. Dispersion Throughout the Zone 

The calculated phonon dispersion curves are displayed 
in Fig. 1. The only noticeable softening is seen midway 
along the F-L direction (transverse T branches), for the < 
110 > polarized 71 branch along the F-K direction empha- 
sized previously, [8] and possibly softening in the T branches 
around the X point. We had identified previously that the 
Q = (|, |, 0)^ 7i phonon (near K) causes large band shifts 
with atomic displacement {Ssk/SuQ « 5 eV/A) near the 
FS necks, reflecting strong EP coupling. 

The corresponding 'mode A' values Aq_^ are presented 
in Fig. 2. The Xq^i, values for the 71 branch near K, and 
also near the midpoint of the F-L line are larger compared 
to other areas in the BZ, as expected from the phonon 
kinks. In addition, we notice that the Xq^i, values for all 
the three branches are large all along F-X; this coupling 
is not immediately obvious from the phonon dispersion in 
Fig. 1 because it depends only weakly on Q. 

In Fig. 3, W6 displciy the corresponding linewidths ^qz^, 
which are given by [9] 

= 27:coQ.J2\^ki+Q\'^i^k)S{sk+Q). (1) 
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Fig. 2. (color online) Calculated mode Aq^ values, following the 
notation of Fig. 1. The landscape is dominated by peaks in the 
transverse branches near K, a longitudinal peak along F-L, and an 
increase around Q~zero where phase space is limited. Here Aq,^ 
is normalized such that the total A requires a sum over the three 
branches (rather than an average). 

While includes the same Q-specific Fermi surface aver- 
age of the EP matrix elements with respect to the available 
phase space for scattering through wavevector Q as does 
Aq^ , it is proportional to ujqi, rather than inversely propor- 
tional. As a result, "fQi, is a much better indicator of the 
importance of the coupling for Tc. Coupling at high fre- 
quency is much more important than at low frequency for 
Tc, whereas the opposite is true for A. Thus it is not sur- 
prising that the linewidth 'dispersion' in Fig. 3 provides a 
different viewpoint than does the Xqi, dispersion. There is 
still the large contribution from the region around K where 
there is an incipient instability, but it is very limited in Q 
and very sharply structured, not being resolved even by the 
Q mesh we have used. The really impressive region however 
is the C branch along (100) directions, which dominates the 
landscape (at least along symmetry lines). There are also 
important contributions from the C branch along the F-K 
lines and from the T branches along (100). Even the small 
peak in the T branches midway between F and L (where 
the kink in the T branch occurs) is small compared to the 
contribution from the C branch along the same line. Tse et 
al. present 7q at V/Vq = 0.45 GPa (« 50 GPa) which have 
some features [12] in common with Fig. 3, but cannot be 
compared in detail because of the different volumes. 

4. Spectral Functions and Tc 

The phonon frequency distribution F{lj) is shown in Fig. 
4 where it is clear that the coarser Q mesh had given a good 
picture of the density of states. This insensitivity does not 
persist when coupling is included. The Eliashberg spectral 
function a^Fiu) from both meshes is plotted in Fig. 5. For 
this 20 GPa case a large broad peak dominates the 15-35 
meV range, with less intense coupling extending on up to 
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Fig. 3. (color online) Calculated linewidths 7q^, following the nota- 
tion of Fig. 1. The landscape here is dominated by all modes along 
F-X, and the region near K. 
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Fig. 4. (color online) Phonon density of states calculated from the 
12^ Q grid (72 points) and the 24^ grid (413 points), illustrating 
that the representation of the phonon spectrum changes little. The 
substantial changes in 0^(0)) shown in Fig. 5 are discussed in the 
text. 

70 nieV. Compared to its value from the coarser mesh, it is 
much smaller in the 7-20 meV range. The intensity in the 
low frequency region is clearly from the transverse branch 
phonons, while the less impressive (but still substantial) in- 
tensity in the high frequency region is from the longitudinal 
phonon branch. 

The frequency resolved coupling strength a'^(aj), also 
shown in Fig. 5 is dramatically changed by using the finer 
Q mesh. There is now a broad peak around 12 meV, still 
indicative of comparatively stronger coupling in the low 
frequency regime. However, the overriding peak in the 7-16 
meV range has disappeared, showing that is was an artifact 
of the sampling of the volume near the point K using only a 
12"^ mesh. This feature is reminiscent of results calculated 
by Yin et al., [14] who have found an even narrower peak in 
a^(w) growing at lower frequencies (2-5 meV) in elemental 
Y under pressure. This peak also was traced to incipient 
instability for a transverse branch near the zone boundary. 

We now concentrate on the difference in Tc, and the in- 
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Fig. 5. (color online) Comparison of the Eliashberg spectral function 
a^Fi^u)) (inset) and the frequency-resolved coupling strength a'^iu)) 
for fee Li at 20 GPa, for both Q point meshes. 

Table 1 

From the calculated cPF(lo) at 20 GPa for both Q meshes, this table 
provides: the logarithmic, first, and second moments of the frequency; 
the value of A; Tc (K) calculated using /i*=0.13 and /i*=0.20; and 
two simple measures of a^F (see text). 
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terpretation of the distribution of the coupling, that we 
obtain from our denser Q mesh of data (we note that the 
quantities we discuss are still not fully converged). The re- 
sults of this much better zone sampling are displayed in 
Table 1. Startlingly, the value of A has dropped from 3.1 to 
1.0, (which was also obtained by Profeta et al. [7] )reflect- 
ing the fact that the low energy region of a^F was (and 
still is) reliant on large contributions (from near K and per- 
haps elsewhere) that vary sharply with Q. For the deter- 
mination of Tc, the strong reduction in A is considerably 
compensated by increases in the frequency moments. The 
logarithmic moment LOiog increases by a factor of 2.7, and 
the second moment by 40%. At a similar volume, Shi and 
Papaconstantopoulos [13] obtain A = 0.5 from rigid mufRn 
tin calculations. One source of difference is the value {lS) « 
325 K (scaled from the P = Debye frequency which is our 
calculated value of wi. 

We have obtained Tc using the Allen-Dynes equation[15] 



Tc 



1.2 



AiA2ea;p{- 



1.04(H-A) 
A-A<*(1 + 0.62A) 



} 



(2) 



where Ai, A2 are strong coupling corrections that depend 
on A, and the ratio 1^2/^^/09, which is a measure of the 
shape of a^F. These A factors are important for A=3.1 but 
have little effect for A ~1 or less. For the commonly used 
value of /i*=0.13, our more nearly converged results give 
Tc=13 K vs. the earlier estimate of 20 K, only a 35% drop 
even though A decreases by more than a factor of three. 
For /Lf*=0.20, Tc = 3.4 K. These values should be compared 




Fig. 6. (color online) Isosurface plot of = 0.054 meV for the 
longitudinal branch. The box contains the F point at the center and 
at each corner for this foe structure. The amplitude is high inside 
the jack-shaped region midway between F and X, corresponding to 
the large linewidths evident in Fig. 3. [Unfaithful interpolation at 
the edges of the box account for the extra structure in those regions; 
the box edges are also F-X-F lines.] 

with the nearly hydrostatic value Tc = 6 K at 20 GPa found 
by Deemyad and Schilling. [2] 



5. Analysis and Summary 

This apparent weak dependence of Tc on A reflects the 
observation, made several times in the past, [16] that, for 
providing an estimate of with averages of a^F , the com- 
bination of A and some frequency moment separately is not 
the best choice. This observation is connected with the ob- 
servation that coupling at low frequency, which strongly 
enhances A, has only a weak effect on Tc. Carbotte and 
collaborators[16] argued that the product (with oji be- 
ing the first moment as defined by Allen and Dynes [15]) is 
nearly proportional to Tc for 1.2< A < 2.4; this quantity is 
proportional to the area under a^F and corresponds to a 
zone integral of ^q^j/ujqu. They were not including strong 
coupling corrections to Tc,[15] however, which were un- 
known at the time of the original work. The strong-coupling 
limit [15] gives Tc oc i/Awf; this product A(jj| corresponds 
to a zone integral of ^Qv We find that, for the fine mesh 
compared to the coarse (see the Table): A decreases by 70%; 
Tc drops by 35% (at /i*=0.13); the Carbotte product de- 
creases by 30%; the strong-coupling limit factor drops by 
only 5%. Hence the Carbotte ratio follows the change in Tc 
best, of these choices. 

The behavior of Xqi, oc ^Qv I^'q^ does not faithfully rep- 
resent the relative importance of modes for Tc except in 
the weak-coupling limit. In the strong coupling limit (which 
we do not claim is the case for Li) it is 7Q,y that gives the 
better picture. For another perspective we display in Fig. 
6 an isosurface of 7q^ for the £ modes. As expected from 
Fig. 3, it is dominated by contributions near the center of 



the (100) line. The isocontours form a blunt jack-like shape 
although the true symmetry is only four-fold. Our corre- 
sponding plots for Aqi/ are almost like those of Profeta et 
al. [7] except that we seem to find somewhat larger contri- 
butions near the zone center. 

In this paper we have clarified and supported our earlier 
suggestion [8] that denser zone sampling would be necessary 
to obtain accurate numerical results. We find in addition 
that careful sampling of the fine structure in the coupling is 
required even for obtaining a correct understanding of the 
microscopic origin of the pressure- induced increase in Tc. In 
particular, we call into question the view (easily surmised 
from the earlier work[7,8] even if not so claimed) that the 
impressive value of Tc in Li under pressure arises either 
from soft phonons near K, or from a box- like set of regions 
in momenta along (?7,?y, C) with 77 ^ 0.7. These questions, 
and the strong polarization dependence of matrix elements, 
will be addressed in further work. 
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